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CN \ Abstract 

■ We consider the fractional generalizations of Liouville equation. The normalization 
pLn '. condition, phase volume, and average values are generalized for fractional case. The 
QQ ' interpretation of fractional analog of phase space as a space with fractal dimension and as 
CN ■ a space with fractional measure are discussed. The fractional analogs of the Hamiltonian 
^ systems are considered as a special class of non-Hamiltonian systems. The fractional 
Q \ generalization of the reduced distribution functions are suggested. The fractional analogs 

' of the BBGKI equations are derived from the fractional Liouville equation. 

d ■ 

■ 1 Introduction 

^ Fractional integrals and derivatives have many applications in statistical mechanics and kinetics 

^ [1]. G.M. Zaslavsky [2, 1, 3] proved that the chaotic dynamics can be described by using the 

\C> Fokker-Planck equation with the coordinate fractional derivatives. The fractional Zaslavsky's 

^ . equation can be used to describe dynamics on fractal [4, 5]. It is known that Fokker-Planck 

O ] equation can be derived from the Liouville and BBGKI equations [6, 7, 8]. The Liouville 

^ \ equation is derived from the normalization condition and the Hamilton equations [9]. In the 

' Hamilton equations we have only the time derivatives. Therefore the usual normalization con- 

• S dition leads to the usual Liouville equation. The BBGKI equations can be derived from the 

^ , Liouville equation and the definition of average value [17, 18, 19, 20]. Therefore Zaslavsky's 

>- ! equation [2, 1, 3] can be derived from a fractional generalization of the Liouville and BBGKI 

k> ' equations. To derive the fractional generalization of the Liouville and BBGKI equations, we 

■ use a fractional generalization of normalization condition [10] and fractional generalization of 
■ — ' the definition of average values. In these generalizations, we use the integrals of fractional 

order. To use these integrals, we must have the physical interpretation of the fractional order 
of integrals. We can consider the fractional integrals as integrals for the function on a frac- 
tional space. In order to use this interpretation of the fractional integrals, we must define a 
fractional space. The interpretation of the fractional space is connected with fractal dimen- 
sion space. We can replace the distribution on fractal with fractal mass dimension by some 
continuous distribution that is described by fractional integrals. This procedure is a fractional 
generalization of Christensen approach [11] that is averaging procedure over fractals. Suggested 
procedure leads to the fractional integration and differentiation to describe the distributions on 
fractal. The fractional integrals allow us to take into account the fractality of the distribution. 
The functions on fractals can be averaged and the distribution on fractal can be replaced by 
some "fractional" continuous distribution. In order to describe the averaged distribution on 
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fractal with non-integer mass dimension, we must use the fractional integrals. Smoothing of 
the microscopic characteristics over the physically infinitesimal volume transforms the initial 
distribution on fractal into distribution on space with measure that uses the fractional inte- 
grals. The order of fractional integral is equal to the fractal mass dimension. The consistent 
approach to describe the distribution on fractal is connected with the mathematical definition 
the integrals on fractals. In Ref. [12], was proved that integrals on net of fractals can be 
approximated by fractional integrals. In Ref. [10], we prove that fractional integrals can be 
considered as integrals over the space with fractional dimension up to numerical factor. We use 
the well-known formulas of dimensional regularizations [16]. Note that almost all systems with 
fractional phase space are non-Hamiltonian dissipative systems in the usual phase space {q,p). 
Therefore we have the other interpretation of the fractional phase space. This interpretation 
follows from the fractional measure of phase space [10] that is used in the fractional integrals. 
The fractional phase space can be considered as phase space that is described by the fractional 
powers of coordinates and momenta. Using this phase space we can consider wide class of 
non-Hamiltonian systems as generalized Hamiltonian systems. In this case, the fractional nor- 
malization condition and the fractional average values are considered as condition and values 
for the generalized Hamiltonian systems that are non-Hamiltonian systems in the usual phase 
space. 

In Sec. 2, we consider the fractional generalization of normalization condition. In Sec. 3, 
we derive the fractional generalization of continuity equation and Liouville equation. In Sec. 
4, we consider the physical interpretations of fractional phase space as a space with fractal 
dimension. In Sec. 5, the fractional phase space is considered as space with fractional measure. 
The fractional generalization of the Hamiltonian systems are suggested. In Sec. 6, we consider 
the fractional generalization of average values. In Sec. 7, the fractional generalization of BBGKI 
equations are derived from the fractional average value and the fractional Liouville equation. 
Finally, a short conclusion is given in Sec. 8. 



2 Fractional Generalization of Normalization Condition 

Let us consider a distribution function p{x,t) for x in 1-dimensional Euclidean space E^. Let 
p{x,t) e Li{E^), where t is a parameter. Normalization condition has the form 

/+0O 
p{x, t)dx = 1. 
-oo 

This condition can be rewritten in the form 

/y r+oo 
p(x,t)dx+ / p(x,t)dx — l, (1) 
-oo Jy 



where y G (— oo, +oo). 

Let p{x,t) G Lp{E^), where 1 < p < 1/a. Fractional integrations on {—oo,y) and {y,+oo) 
are defined [14] by the equations 

, 1 fv p(x.t)dx /^n, ^/ X 1 f~^°° p(x,t)dx 
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Using these notations, we rewrite Eq. (1) in an equivalent form 

{Ilp){y,t) + {l\p){y,t)^l. 

The fractional analog of normalization condition (1) can be represented by the following equa- 
tion 

Equations (2) can be rewritten in the form 

{I^p){y,t) = — - / x^-'p{yTx,t)dx. (3) 
i [a) Jo 

This leads to the fractional generalization of the normalization condition 

p{x,t)dpa{x) = 1, (4) 



where we use the following notations 

p{x,t)^T^p^ -(p{y-x,t)+p{y + x,t)j, dpa{x) ^ ^^^^ dx. (5) 

The normalization in the phase space is derived by analogy with a normalization in the 
configuration space. The fractional normalization condition in the phase space is 

/+00 r+oo 
/ p{Q^P,t)dpa{q,p)^'i-, (6) 
-oo J —oo 

where dpa{(l,p) has the form 

dpa{q,p) = dpaiq) A dpa{q) = T^n. V dq A dp. (7) 

The distribution function p{q, p, t) in the phase space is defined by 

p{q,p,t)^TgTpp{q,p,t), (8) 
where the operators Tg and Tp are defined by the equation 

TxJ{-;Xk,--) = ^[f{...,x'^-Xk,...) + f{...,x'^ + Xk,...)). (9) 
This operator allows us to rewrite the distribution functions 

p{q,P,t) = ^{p{q' - q,p' - P,t) + p{q' + q,p' - P,t) + p{q' - q,p' + P,t) + p{q' + q,p' + p,t 
in the simple form p{q,p,t) = TqTpp{q,p,t). 



3 



3 Fractional Generalization of Liouville equation 

3.1 Fractional Generalization Continuity Equation 

Let us consider a region Wq for the time t — 0. In the Hamilton picture, we have 

JWt JWo 

Using the replacement of variables Xt = Xt{xo), where xq is a Lagrangian variable, we get 

-idxt 



p{xt,t)\xt\" ^^dxo = [ p(a;o,0)|a;o|" ^dxo. 

OXq J Wo 



I Wo 

Since Wq is an arbitrary region, we have 

p{xf, t)dna{xt) = p{xo, 0)diJ,a{xo). (10) 
This condition leads us to the equation 

p{x„t)\x,r-'^^^ p{xo,o)\xor\ 

Differentiating this equation in time t, we obtain 

As the result we have the fractional generalization of the continuity equation 

^Ei^ + n^^:,^^t)p{xt,t)^o, (11) 

where the omega function 

describes the velocity of phase volume change. Eq. (11) is a fractional continuity equation in 
the Hamilton picture. If the equation of motion has the form 

dt 

then the function is defined by 

r, / ^\ ^ (i \ la-i ,7 ^^t\ a-ldxt d dxt {a-l)Ft dFt , . 

dt ^ oxo ' Xt dt oxt dt Xt oxt 
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3.2 Fractional Continuity Equation for Phase Space 

Using phase space analog of Eq. (10) in the form 

PtdHa{quPt) ^ PodlIa{qo,Po), 

we get the relation 

Pt \ ^ ^^'^ = po ° ° , dqo A (ipo- (13) 

Let us use the well-known transformation 

dqt A dpt = {qt,Pt}odqo A dpo, (14) 
where {qt,Pt}o is Jacobian which is defined by the determinant 

r , .diquPt) , , / dqkt/dqio dqu/dpiQ \ 

As the result, we have condition (13) in the form 

Pt\qtPt\'^~^{qt,Pt}o = koPoT'Vo- (15) 
This equation can be rewritten in more simple form 

Pt{q?,Pt}o^a'\qoPor'po. (16) 

We use the following notation for fractional power of coordinates and momenta 

= /3(x)(x)" = s^n(x)|x|", (17) 

where j3{x) = {sgn{x))°''^^ . The function sgn{x) is equal to +1 for x > 0, and —1 for x < 0. The 
total time derivatives of Eq. (16) lead us to the fractional generalization of Liouville equation 
in the form 

^ + ^^aP = 0, (18) 

where the omega function VL^ is defined by 

= {qt,P^]o'j^{qt.P't}o = |M5r,P?}o- (19) 



In the usual notations, we have 



= ^In det^^^. (20) 
dt d{qo,po) 



Using well-known relation In det A = Sp In A, we can write the omega function fla in th^ form 



5 



where { , is the fractional generahzation of the Poisson brackets in the form 

dA dB dA dB 



In the general case (a 7^ 1), the function Qq, is not equal to zero (O^ 7^ 0) for the systems 
that are Hamiltonian systems in the usual phase space. If a = 1, we have Qq, 7^ only for 
non-Hamiltonian systems. If the Hamilton equations have the form 

^ = ^-^ = n<iuPt\ (21) 

then the omega function is defined by 

= (a - l)(q-^K{q,p) +p-'F{q,p)) + {K,p}, + {q,F}^. (22) 

This relation allows to derive for all dynamical systems (21). It is easy to see that the usual 
nondissipative system 

dq P dp 

has the omega function 

^a{q:P) = {a - l){mqp)~^{p^ + mqf{q)), 
and can be considered as a dissipative system in the fractional phase space {q°',p°')- 

3.3 Fractional Liouville Equation for N-particle System 

Let us consider A'"-particle system. Suppose /c-particle is described by the generalized coordi- 
nates qfc = (gfci, qkn) ^-nd generalized momenta = {pki, ■■■,Pkn) that satisfy the Hamilton 
equations in the form 

^ = K,(q",p"), ^ = F,(q",p",t). (24) 

Here, we use the notation (17) for fractional power of coordinates and momenta. The evolution 
of N-particle distribution function pjv is described by the Liouville equation. Using the fractional 
normalization condition 

r[l,...,7V]p^(q,p,t) = l, (25) 

we can derive the fractional Liouville equation [10] for N-particle distribution function in the 
form 

^ + Q^p^ = 0, (26) 

where d/dt is a total time derivative 

d d dqk d ^ d,pk d 



dt dt dt (9qfc dt dpk 
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Here, we use the following notations 



a=l 

Using Eq. (24), the total time derivative can be written for the fractional powers in the form 

d d ^ d ^ d 
dt dt aq^ 

The omega function is defined by the equation 

TV n 

= EEii^aMc + Uta. }a). (28) 

k=l a=l 

Here, we use the following notations for the brackets 

_^^fdA_dB__dA_dB_^ 

fc=l a=l "yfca "Pka "Pka 

Using Eqs. (26), (28) and (27), we can rewrite the Liouville equation in the equivalent form 

^ = CnPn, (30) 
where >Cjv is Liouville operator that is defined by the equation 

^"""""S^^^+^sr)- '''' 



4 Fractional Space as Space with Fractal Dimension 

The interpretation of the fractional space can be connected with fractal mass dimension. Fractal 
dimension can be best calculated by box counting method which means drawing a box of size 
R and counting the mass inside. The mass fractal dimension [4, 5] can be easy measured. The 
properties of the fractal like mass obeys a power law relation 

M{R) = kR'^, (32) 

where M is the mass of fractal, i? is a box size (or a sphere radius), and a is a mass fractal 
dimension. Amount of mass inside a box of size R has a power law relation (32). 

Let us consider the region W in 3-dimensional Euclidean space E^. The volume of the 
region W is denoted by y(l^). The mass of the region W of the distribution on fractal is 
denoted by M{W). The fractality means that the mass of this distribution in any region W of 
Euchdean space increase more slowly that the volume of this region E^. For the ball region 
of the distribution on fractal, this property can be described by the power law (32), where R 
is the radius of the ball W. 
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The distribution on fractal is called homogeneous if the power law (32) does not depends 
on the translation and rotation of the region W. The homogeneity property of the distribution 
can be formulated in the form: For all regions Wi and W2 such that the volumes are equal 
V^(H^i) = K(W^2)) we have that the mass of these regions are equal M{Wi) = M{W2). In order 
to describe the homogeneous distribution on fractals, we must use the continuous distribution 
such that fractality and homogeneity properties can be realized in the form: 

(1) Fractality: The mass of the ball region W for the distribution on fractal obeys a power law 
relation (32). In the general case, we have the scaling law relation 

where XW = {Xx, x e W}. 



(2) Homogeneity: The local density of homogeneous fractal distribution is translation and 
rotation invariant value that have the form p{x) — Po — const. 

We can realize these requirements by the fractional generalization of the equation 

M3(W) = f p{x)d^x. (33) 



IW 

Let us define the fractional integral in Euclidean space in the Riesz form [14] by the equation 



{rp){y) = / pix)dpaix), (34) 



w 



where d/iaix) — Cs{a,x,y)d^x, and 



, , 23-"r(3/2), ,^_3 , , 
C3[a,x,y)^ \x-y\ , \x - y\ ^ 
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r(«/2) \ ^1 

The point y & W \s the initial point of the fractional integral. We will use the initial points in 
the integrals are set to zero {y = Q). The numerical factor in Eq. (34) has this form in order to 
derive usual integral in the limit a ^ (3 — 0). Note that the usual numerical factor 73^^(0;) = 
r(l/2)/2"7r3/2r(a/2), which is used in Ref. [14], leads us to 73-^(3 - 0) = r(l/2)/2V/2r(3/2) 
in the limit a ^ (3 — 0). 

Using notations (34), we can rewrite Eq. (33) in the form M^iyV) — {Pp){y). Therefore 
the fractional generalization of this equation can be defined in the form 

M^{W) = {rp){y) = )' ^ / p{x)\x - yr'd'x. (35) 

1 {a/ 2) Jw 

If we consider the homogeneous fractal distribution {p{x) = po = const) and the ball region 
W — {x : \x\ < R}, then we have 



r(Qi/2) Jw 
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Using the spherical coordinates, we get 



Ma{W) 



7r25-"r(3/2) 

r(«/2) 



25-"7rr(3/2) 
ar(a/2) 



As the result, wc have M(Vr) ~ i?", i.e., we derive Eq. (32) up to the numerical factor. 
Therefore the distribution on fractal with non-integer mass dimension a can be described by 
fractional integral of order a. 

Note that the interpretation of the fractional integration is connected with fractional di- 
mension [10]. This interpretation follows from the well-known formulas for dimensional regu- 
larizations [16]: 



up to the numerical factor r{a/2)/{2n°'^'^r{a)). 

5 Fractional Space as Space with Fractional Measure 

The interpretation of the fractional space is connected with the fractional measure that is used 
in the fractional integrals. The parameter a defines the space with the fractional measure 
(volume) of the region W. It is easy to prove that the velocity of the fractional measure 
(volume) change is defined by the omega function (12). 

5.1 Fractional Phase Volume for Configuration Space 

The usual phase volume of the region W — {x : x E [a;b]} in Euchdean space is defined by 




(36) 



Using Eq. (36), we get that the fractional integral 




can be considered as integral in the fractional dimension space 




(37) 




(38) 



where y e [a; b]. Using the fractional integrals [14] in the form 




we get the phase volume (38) in the equivalent form 



MW)^{lUmy) + {lll){y). 



(39) 
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The fractional generalization of the phase volume can be defined by 

lJ^a{W) = {i:^l){y)-r{It^){y). (40) 
The fractional phase volume integral can be represented by the equation 

^J^a{W) = / g{a)d^ia{x). (41) 

J-ip-y) 

Here, we use the notations 

= ^(5^ = M") = 5(1 + (fr^) ). (42) 

We use the notation (17) for fractional power of coordinates. 
5.2 Fractional Phase Volume for Phase Space 

The fractional measure for the region W of 2n-dimensional phase space can be defined by the 
equation 

l^aiW) = / g{a)dna{q,p), (43) 
Jw 

where diJ,a{q,p) is a phase volume element 

and g{a) is a numerical multiplier. If the region W of the phase space is defined by G 
and pk e E'^, then g{a) = 1/4''. If this region is defined by Qk G [qakiqbk] and Pk G [paklPbk], 
then 

It is easy to see that the fractional measure depends on the fractional powers of coordinates 
and momenta (5",]?"). 

For example, the phase volume for the two-dimensional phase space has the form 

J-{qb-y')J-(s>b-y') {aV{a)Y J-{q,-y') J-(j>,-y') (r(Q;))2 

Note that the volume element of fractional phase space can be realized by fractional exterior 
derivatives [13] that are defined by 

n Qa n Qa 

= S "^'hdiqk-vk))- ^ S '^''hd{Pk-y'k)r 

For example, the two-dimensional phase space is defined by 
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5.3 Fractional Phase Volume Change 



The interpretation of the fractional phase space is connected with the fractional measure of 
phase space. The parameter a defines the space with the fractional phase measure (43) and 
(44) . It is easy to prove that the velocity of the fractional phase volume change is defined by 
the equation 



dt Jw 
where the omega function Qq. has the form 



Here, we use the fractional brackets: 

where we use notations (17). 

The form of the omega function allows us to consider new class of the systems that are 
described by the fractional powers of coordinates and momenta. The system can be called a 
fractional dissipativc system if a fractional phase volume changes, i.e., fi„ ^ 0. The system 
which is a nondissipative system in the usual phase space, can be a dissipative system in the 
fractional phase space. The usual nondissipative systems (23) are dissipative in the fractional 
phase space. 



5.4 Fractional Generalization of Hamiltonian Systems 

Fractional phase space can be considered as a phase space of the systems that are described by 
the fractional powers of coordinates and momenta. Let us consider the class of non-Hamiltonian 
systems that are described by the fractional powers of coordinates and momenta. A system is 
called a fractional system if the phase space of the system is described by the fractional powers 
of coordinates and momenta 

qt = m{qkT = sgn{qk)\qk\'', pt = m{PkT = sgn{pk)\pk\" ■ (49) 

Here A: = 1, ...,n, and I3{x) is defined by Eq. (17). We can consider the fractional systems in 
the usual phase space (g,p) and in the fractional phase space (q'",p"). In the second case, the 
equations of motion for the fractional systems have more simple form. Therefore we use the 
fractional phase space that is a space with the fractional measure that is used in the fractional 
integrals. We consider the fractional power of the coordinates as a convenient way to describe 
systems in the space with measure that is defined by fractional integrals. 

A classical system (in the usual phase space) is called Hamiltonian if the right-hand sides 
of the equations 

% = i^.(5,p), '^ = F,{q,p) (50) 
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satisfy the following Helmholtz conditions [15]: 

dKk dKi dKk dFi dF^ dFi , , 

dp; dpfc dqi dpk oqi oqk 

In this case, we can rewrite Eqs. (50) in the form 

dqk dH dpk dH 

dt dpk dt dqk 

Using the Poisson brackets { , }i, we can rewrite these equations in an equivalent form 



The fractional generalization of the Hamiltonian system is described by the equation 

dq^ _ dH dp^ __dH 
~dt~'d^' 'df~~dq^' 



(52) 



where if is a function that can be considered as a fractional analog of the Hamiltonian. Note 
that the function H such that dH/dt = is the invariant of the motion. Using the brackets 
(48), we can rewrite Eq. (52) in the equivalent form 

These equations describe the system in the fractional phase space {q°',p"). For the usual phase 
space {q,p), the fractional Hamiltonian systems (52) are described by the equations 

dqk ^ (gfcPfc)^"" dH dpk ^ (gfcPfc)^"" dH 
dt o? dpk dt oP- dqk 

The fractional Hamiltonian systems (54) are non- Hamiltonian systems in the usual phase space 
{q,p)- It is easy to prove that the Helmholtz conditions (51) are not satisfied. Therefore 
fractional Hamiltonian system (54) is a non-Hamiltonian system in the usual phase space {q,p). 
The fractional phase space allows us to write the equations of motion for the non-Hamiltonian 
systems (54) in the generalized Hamiltonian form (52). 

The omega function for the system (50) in the usual phase space {q,p) is defined by the 
equation 

If the omega function is negative i7 < 0, then the system is called a dissipative system. If 7^ 0, 
then the system is a generalized dissipative system. For the fractional Hamiltonian systems 
(54), the omega function (55) is not equal to zero. Therefore the fractional Hamiltonian systems 
are the general dissipative systems in the usual phase space. 

The function H = H{q'^,p'^) can be considered as a fractional analog of the Hamiltonian 
function. For example, we can use 

n 2a 

H{q-,p-) = J2^ + U{q'^). (56) 
fe=i 
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It is easy to see that fractional systems with Hamiltonian (56) lead us to the non-Gaussian 
statistics. The interest in and relevance of fractional kinetic equations is a natural consequence 
of the realization of the importance of non-Gaussian statistics of many dynamical systems. 
There is already a substantial literature studying such equations in one or more space dimen- 
sions. 

Note that the classical dissipative non- Hamiltonian systems can have canonical Gibbs distri- 
bution as a solution of stationary Liouville equations for this dissipative system [23] . Using the 
methods [23], it is easy to prove that some of fractional dissipative systems can have fractional 
generalization of the canonical Gibbs distribution in the form 

p{q,p) = Z{T)exp — , 

as a solution of the fractional Liouville equations 

Here the function H{q",p°') is defined by (56). 

6 Fractional Generalization of Average Values 

The usual average value for the configuration space 

/+00 
f{x,t)p{x,t)dx (58) 
-oo 

can be written in the form 

ry r+oo 
< f >i= f{x,t)p{x,t)dx + f{x,t)p{x,t)dx, (59) 



where y e (— oo,-|-oo). Using the notations (2), we can rewrite the average value (59) in the 
form 

< / >i {y,t) = {Ilfp){y,t) + {Iifp){y,t). 
The fractional generalization of this equation is defined by 

< / >« {y,t) = {nfp){y,t) + {rfp){y,t). (60) 

We can rewrite Eq. (60) in the form 

/ + C10 
T,{fp){x,t)dpa{x), (61) 
-oo 

where we use 

d.„W = ^ = ^. r./ = !(/(, + + (62) 
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and is defined by Eq. (17). 

We can define the integral operator 7" by the equation 

/ + 00 
f{x)dfia{x). (63) 
-oo 

In this case, the fractional generalization of average value (61) can be written in the form 

<f>a {y,t)^i^TJ{x)p{x,t). 

Let us consider /c-particle that is described by generalized coordinates qfc = (g^i, qkn) and 
generalized momenta = (pki, ■■■iPkn)- We can define the phase space integral operator for 
/c-particle by the equation 

J<^\U\ — f" f" fa f" 

For the A^-particle system, we use the operators 

7"[1, A^] = r[l, N\ = T\V\...T\N\ 

where the operator T[A;] is defined by the relation 

Here the operator T^u is defined by Eq. (62). 

Using these notations, we can define fractional analog of the average values for the phase 
space of A^-particle system by the relation 

< / (y, t) = /"[I, N]T[1, iV]/p^. (64) 

In the simple case {N — m — 1), the fractional average value (64) is defined by the equation 

/'+00 r+co 
/ dfia{q,p) TqTpf{q,p)p{q,p,t). (65) 
-oo J —oo 

Note that the fractional normalization condition is a special case of this definition of the average 
value < 1 >a (y, t) — 1. 

7 Fractional Generalization of BBGKI Equations 

The state of N-particle system is described by N-particle distribution function 

PAr(q, P, t) = p(qi, pi, Qat, Pat, t). 
We use the tilde distribution functions 

PAr(q, p, t) = r[l, N]pN{q, p, t). (66) 

The fractional generalization of 1-particle reduced distribution function pi can be defined by 
the equation 

pi(q,p,i) =p(qi,Pi,i) =/"[2,...,7V]p^(q,p,i). (67) 
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Obviously, that 1-particle distribution function satisfies the fractional normalization condition 

r[l]pi(q,p,t) = l. (68) 

The BBGKI equations [17, 18, 19, 20] are equations for the reduced distribution functions. 
These equations can be derived from the Liouville equation. Let us derive the fractional gen- 
eralization of the first BBGKI equation from the fractional Liouville equation (30). 

In order to derive the equation for 1-particle distribution function pi we differentiate Eq. 
(67) with respect to time: 

Using the Liouville equation (30) for A'"-particle distribution function p^, we have 

^^n2,...,N]C^Mq,p,t). (69) 



Substituting (31) in Eq. (69), we get 

^Pi faro An ^i'^(I^fePAf) , 9{FkpN) 
k=l 



Let us consider in Eq. (70) the integration over and for k-particle term. Since the 
coordinates and momenta are independent variables, we can derive 

rM4(K*P.) = ^(K.^„):>0. (71) 

Here, we use that the distribution p^r in the limit ±00 is equal to zero. It follows from 

the normalization condition. If the limit is not equal to zero, then the integration over phase 
space is equal to infinity. Similarly, we have 

d 



^dpt'^^'-^J aria) 

Then all terms in Eq. (70) with k = 2, ...,n are equal to zero. We have only term for k = 1. 
Therefore Eq. (70) has the form 

f ^-M2....,iVl(^.^). (72) 

Since the variable qi is an independent of q2,...,qiv and p2,...,PAr, the first term in Eq. 
(72) can be written in the form 

L ' ' J 9qf ^ ' ' aqf 

The force Fi acts on the first particle. This force is a sum of the internal forces Fi^ = 
F(qi, Pi, qfc, Pfc, t), and the external force F^ = F'^(qi, pi, t). In the case of binary interactions, 
we have 

N 

F,=Fl + J2F,k. (73) 

k=2 
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Using Eq. (73), the second term in Eq. (72) can be rewritten in the form 

We assume that distribution function is invariant under the permutations of identical particles 
[21]: 

PN{.:,qk,Pk,--,(li,Pi,--,t) = pjv(--,qi,Pi,-,qik,Pfe, •••.^)- 

In this case, the N-particle distribution function is a symmetric function for the identical 
particles and we have that all {N — 1) terms of sum (74) are identical. Therefore the sum can 
be replaced by one term with the multiplier (A^ — 1): 

E/12,...,A^] ^(f,,p^) = (7V-1)7"[2,...,7V] ^(f^^Pa.). (75) 
Using 7°^[2, A^] = 7°^ [2] 7° [3, A^], we can rewrite the right hand side of Eq. (75) in the form 

Here, we use definition of 2-particle distribution function p2- This distribution is defined by the 
fractional integration of the A'"-particle distribution function over all qfc and Pk, except k — 1,2: 

P2 = p(qi, Pi, q2, P2, t) = r[3, N]pN{q, p, t). (77) 

Since pi is independent of q2, P2, we can change (76) the order of the integrations and the 
differentiations: 

Finally, we obtain the equation for one-particle reduced distribution function 

dpi , .9(Kipi) , djFlpi) 

-dr^-d^^-d^ = ^^p'^- ^^^^ 

Here /(P2) is a term with 2-particle reduced distribution function 

7(p2) = -(Ar-l)^r[2]Fi2P2. (79) 

Therefore the fractional generalization of the first BBGKI equation has the form 

dpi 



dt 

where £1 is 1-particle Liouville operator 



+ /(pa), 



^'P' - "^^f 
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The physical meaning of the term /(P2) is following. The term 7(p2)(i/Xa(q, p) is a velocity of 
particle number change in 4m-dimensional elementary phase volume o?//a(qi, P2, ^2, P2)- This 
change is caused by the interactions between particles, li a — 1, then we have the first BBGKI 
equation for non-Hamiltonian systems. 

Let us consider the particles as statistical independent systems. In this case, we have 

P2(qi, Pi, q2, P2, t) = pi(qi, pi, i)pi(q2, P2, t). (81) 
Substituting (81) in (79), we get 

I(P2) = -^pir[2]Fi2Pi(q2,P2,i), 

where pi — pi(qi, Pi,t). As the result, we have the effective force 

F^-^^(qi,Pi,t) = r[2]Fi2Pi(q2,P2,t). 
In this case, we can rewrite the term (79) in the form 

Substituting (82) in Eq. (78), wc get 

This equation is a closed equation for 1-particle distribution function with the external force Ff 
and the effective force F'^-'^-'^. Equation (83) is a fractional generalization of the Vlasov equation. 

Let us differentiate this equation (77) that defines the two-particle reduced distribution 
function p2. The fractional Liouville equation allows us to derive equation for 2-particle reduced 
distribution function p2 in the form 

^= AP2 + >C2P2+A2P2 + /(P3), (84) 

ot 

where Ci and £2 are 1-particle Liouville operators (80) and C12 is 2-particle Liouville operator 
that is defined by equation 

d d 

A2P2 = ^(Fi2P2) + ^(F2iP2), 



and /(pa) is a term with the 3-particle reduced distribution 

,(,3). ("-^("-^) rl3l(g(g^.g(ggM). (85) 



The derivation of Eq. (84) is the analogous to the derivation of Eq. (78). It is easy to see that 
the system of Eqs. (78) and (84) are not closed. The system of these equations for the reduced 
distribution functions can be called the fractional generalization of the BBGKI equations. 
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8 Conclusion 



In this paper, we consider the fractional generahzations of Liouvillc and BBGKI equations. 
The normahzation condition, phase volume, and average values arc generalized for fractional 
case. These generalizations lead us to the fractional analog of phase space. The space can be 
considered as a fractal dimensional space. The physical interpretation of the fractional phase 
space is discussed. The fractional generalization of average values is derived. 

In this paper the fractional analogs of the BBGKI equations arc derived. In order to derive 
these analogs we use the fractional Liouville equation [10], we define the fractional average values 
and the fractional reduced distribution functions. The fractional analog of Vlasov equation is 
considered. 

The fractional Liouville, BBGKI and Vlasov equations are better approximation than its 
classical analogs for the systems with the fractional phase spaces. For example, the systems 
that live on some fractals can be described by the suggested fractional equations. Note that 
the fractional system is non-Hamiltonian systems. The fractional harmonic oscillator is an 
oscillator in the fractional phase space that can be considered as a fractal medium. Therefore 
the fractional oscillator can be interpreted as an elementary excitation of some fractal medium 
with non-integer mass dimension. 

It is not hard to prove that the fractional systems are connected with the non-Gaussian 
statistics. That the dissipative and non-Hamiltonian systems can have stationary states of 
the Hamiltonian systems [22]. Classical dissipative and non-Hamiltonian systems can have 
the canonical Gibbs distribution as a solution of the stationary Liouville equations for this 
dissipative system [23]. Using the methods [23], it is easy to prove that some fractional dis- 
sipative systems can have fractional analog of the Gibbs distribution (non-Gaussian statistic) 
as a solution of the fractional Liouville equations. Using the methods [23], it is easy to find 
the stationary solutions of the fractional BBGKI equations for the fractional systems. Note 
that the interest in and relevance of fractional kinetic equations is a natural consequence of the 
realization of the importance of non-Gaussian statistics of many dynamical systems. There is 
already a substantial literature studying such equations in one or more space dimensions. 

Note that the quantization of the fractional systems is a quantization of non-Hamiltonian 
dissipative systems. Using the method, which is suggested in Refs. [24, 25], we can realize the 
Weyl quantization for the fractional systems. The suggested fractional Hamilton and Liouville 
equations allow us to derive the fractional generalization for the quantum systems by methods 
suggested in Refs. [24, 25]. 

The fractional BBGKI equations can be used to derive the Enskog transport equations. 
The fractional analog of the hydrodynamic equations can be derived from the first fractional 
BBGKI equation. It is known that the Fokker-Planck equation can be derived from the BBGKI 
equations [6]. The fractional Zaslavsky's equation [1, 2] can be derived from the fractional 
generalization of the BBGKI equation. 
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